Beam shaping by a layered structure with left-handed materials 
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We analyze transmission of a layered photonic structure (a one-dimensional photonic crystal) 
consisting of alternating slabs of two materials with positive and negative refractive index. For the 
periodic structure with zero averaged refractive index, we demonstrate a number of unique properties 
of the beam transmission observed in strong beam modification and reshaping. 
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Recent experimental results confirmed the existence of 
a novel type of composite materials, the so-called left- 
handed (LH) metamaterials (or materials with both neg- 
ative permittivity and permeability), which possess the 
property of negative refraction. Such materials were sug- 
gested theoretically a long time ago [||, but they have 
recently attracted much attention due to their experi- 
mental realization and recent 'hot' debates on the use of 
a LH slab as a perfect lens to focus both propagating 
and evanescent waves The perfect lens concept was 
first suggested by Pendry who demonstrated that a 
slab of a lossless negative-refraction material can provide 
a perfect image of a point source. Although the perfect 
image is a result of an ideal theoretical model employed 
in the analysis of Ref . |^ , the resolution limit was shown 
to be independent on the wavelength of electromagnetic 
wave (but can be determined by other sources such as 
loss, spatial dispersion, etc.), and it can be better than 
the resolution of a conventional lens Q . 

Multilayered structures containing LH materials can 
be considered as a sequence of the perfect lenses and, 
therefore, they are expected to possess novel and unique 
transmission properties. Such multilayered structures 
have been investigated theoretically for calculating the 
transmittance or reflectance in the Bragg regime ^. 
More recently, it was shown that a one-dimensional stack 
of layers with alternating right- and left-handed materi- 
als with zero averaged refractive index displays a narrow 
spectral gap in the transmission ^, ^, which is quite 
different from a Bragg reflection gap. 

In this Letter, we study the properties of layered pho- 
tonic structures consisting of alternating slabs of positive 
and negative refractive index materials (see Fig. |l|), and 
demonstrate unusual angular dependencies for the trans- 
mission of such slabs when the averaged refractive index 
is close to zero. We demonstrate how these properties 
can be employed for strong beam reshaping. 

We consider a one-dimensional photonic crystal formed 
by alternative slabs, as schematically shown in Fig. |^. 
We assume that the periodic structure is made from the 
materials of two types: the slabs of the width b are made 
of a LH material, and they are separated by the right- 
handed (RH) layers of the width a. Such a structure can 
be treated as an array of Pendry's perfect lenses with the 




FIG. 1: Schematic of a multilayered structure consisting of 
slabs with alternating right- and left-handed materials. 



variation of the refractive index in one structural period. 
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where and ni are the positive and negative indices 
of refraction of the RH and LH materials, respectively. 
First, we consider waves propagating in the {x, z) plane 
with the wave vectors k = (kx,0,kz), which are TE- 
polarized, i.e. they have the only component E = Ey 
described by the Helmholtz-type equation. 



AE+^n^z)E- 



1 dfi 
/i(x) dx 



E = 0, 



(2) 



where A is the two-dimensional Laplacian. In an infinite 
periodic structure the propagating waves have the form 
of Block modes, which electric field envelopes satisfy the 
periodicity property, E{z + A) = E{z) exp{iKh). Here 
Kb is the normalized Bloch wave number, and it can be 
found as a solution of the standard equation for two- 
layered periodic structures (see, e.g., Ref. B), 
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Here krz and kiz are the z-components of the wavevec- 
tor in the RH and LH media, respectively, and Pr,i — 

\J <^r,ll [ir,U COs{6r,l), whcrC 0r,l = tan^^ {k(r,l)x/k{r,l)z) 

are the propagation angles in the corresponding media. 
For real Kb the Bloch waves are propagating; complex 
Kb indicate the presence of hand gaps, where the wave 
propagation is inhibited. 
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FIG. 2: Band-gap structure on the parameter plane {Lu,kx) 
with gaps shaded, (a) Transmission band corresponds to a 
normal incidence; (b) Transmission band corresponds to an 
oblique incidence. A dotted line is the frequency uq. Insets 
show the beam transmission coefficients at the frequencies 
marked by dashed lines. 



In comparison with the conventional case of the peri- 
odic structure made of two different RH materials, the 
structure with LH layers shown in Fig. ^ possesses a dis- 
tinctive bandgap of a new type when Jfcirjal = l/cj^^j, the 
property recently revealed in Refs. [[7| Q]- The inherent 
feature of the LH materials is their frequency dispersion, 
so that in most of the cases we can find the characteris- 
tic frequency luq such that rir = |ni|. Then, a periodic 
structure consisting of the alternating LH and RH slabs 
of the same thickness (a = 6) will formally possess a 
bandgap for all angles of incidence, with the only excep- 
tion of the 100% transmission resonances, which appear 
when the arguments of the sinus functions in Eq. (j^) 
are equal to the whole multiple of tt. Such highly un- 
usual transmission can happen because the "zero-index" 
bandgap does not depend on the optical period of the 
structure, whereas usual Bragg-reflection and transmis- 
sion resonances are highly sensitive to variations of the 
period and the incident angle. 

It is important to study the transmission properties 
of a realistic system when the resonant conditions may 
not be exactly satisfied. To be specific, let us consider a 
periodic structure consisting of the layers of LH metama- 
terial separated by air. We assume that the composite 
material possesses the negative refractive index in the mi- 
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FIG. 3: Cross-section intensity profiles of the incident (left), 
reflected (middle), and transmitted (right) beams of various 
shapes: (a) Gaussian, (b-c) Gaussian with the vortex (topo- 
logical) charge four at (b) the normal and (c) oblique incident 
angles. 



crowave region with the effective dielectric permittivity 
and magnetic permeability H, 



e(c.) 



(4) 



where Wp/27r = 10 GHz, tjj27r = 4 GHz, and F = 0.56. 
The frequency range, where both e and /i are negative 
is (4 - 6) GHz. We plot the corresponding band-gap on 
the plane {Lo,kx) [see Fig. ||]. We reveal that for a fixed 
incident angle {kx) there appears a single frequency gap 
inside the transmission region. Quite a different situa- 
tion occurs for the anglular dependence of the transmis- 
sion coefficient at a given frequency, when there exists a 
narrow transmission band inside the bandgap. 

In Fig. |(a), the period is chosen to satisfy the con- 
dition fc^flo = TT (in our case, Oq = 6o ~ 2.88 cm), and 
the resonant transmission of electromagnetic waves of the 
frequency wq occurs at the normal incidence. In contrast. 
Fig. 11(b) corresponds to the case a = = oq + 0.03 cm 
when the transmission resonance occurs for the waves 
at the incidence angle 9: kr cos 6 a — tt. When the pe- 
riod is further increased, multiple transmission bands are 
observed. The width of the transmission band in the k- 
space can be controlled by tuning the frequency in the 
vicinity of wq. The insets in Fig. |^ show the cross-sections 
of the transmission coefficients for the frequencies marked 
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FIG. 4: The Goos-Hanchen shift A of the reflected beam vs 
kx- Solid curves 1,2,3 — numerically calculated shifts for 
the incident Gaussian beams 20, 60, and 100 cm wide, re- 
spectively. Dashed curve — the asymptotic result for wide 
beams. Shaded area corresponds to the transmission band of 
the infinite periodic structure. 

by dashed lines. Thus, one can make a periodic structure 
with the transmission band at the desired angle of inci- 
dence by a proper choice of the system period, while the 
width of the transmission band can be adjusted by tuning 
the frequency. Similarly, the structures with more com- 
plicated transmission properties which include multiple 
rings in the transmission coefficient vs wavevector can be 
obtained. Most importantly, we find that such narrow- 
band angular filters can have almost 100% transparency 
in the transmission region. 

We now use the standard transfer matrix method to 
analyse beam shaping by a finite composite structure 
consisting of 100 layers with the parameters correspond- 
ing to the angular transmission coefficient shown in the 
inset of Fig. ^(b). We consider the transmission of (2-1-1)- 
dimensional beams in the paraxial approximation when 
the angles of incidence are small and polarization effects 
are negligable. For a Gaussian beam incident at a nor- 
mal angle, we observe almost complete reflection if the 
the angular spectrum width is comparable to the width 
of the transmission band. Transmission is only possible 
at the incident angles when the spectra overlap, and un- 
der such conditions the reflected beam has a two-humped 
shape, see Fig. 0(a). On the other hand, vortex beams 



have a ring structure of the angular spectrum and can 
therefore be very effectively transmitted even at the nor- 
mal incidence [Fig. §(b)]. Similar to the case of a Gaus- 
sian beam, we observe that the weak reflected beam has a 
two-hump vortex structure. Small variations of the angle 
of incidence from the normal destroy the structure of the 
reflected and transmitted vortex beams and the results 
depend on the vortex topological charge [see Fig. |3|(c)]. 

A beam reflected by an interface experiences a shift of 
its center relative to the center of the incident beam, and 
this classical effect is known as the Goos-Hanchen shift 
(see, e.g. Ref. ^^)- Reflection from a single LH slab 
results in a negative Goos-Hanchen shift |l^ , and it is in- 
teresting to calculate the beam shift for the layered struc- 
ture. The centers of the incident and reflected beams 
are found as X.^^^ = ^°°^r±_\E,^r\'^ dv^/ \E,^rVdT^, 
where r_L is the transverse coordinate in the plane (x, y). 
For wide beams (i.e. the beams with a narrow spectrum), 
the shift A = — Xj is determined by the gradient of 
the phase ip of the reflection coefficient, A = — Vk^ (j)- 
We calculate the Goos-Hanchen shift for the beam re- 
fiected by the periodic RH/LH structure described above, 
using the parameters corresponding to Fig. ||(b). With 
no loss of generality, we assume that the {x, z) plane is 
the plane of incidence, and therefore the beam is only 
shifted along the x axis. We plot the shifts experienced 
by the Gaussian beam of various widths (see the curves 
1-3 in Fig. ||) . We see that the shift of wide beams indeed 
approaches the the asymptotic result shown with dashed 
line in Fig. |. The shift can be either positive and neg- 
ative, and the absolute value increases dramatically for 
the angles of incidence close to the transmission band. 

In conclusion, we have analyzed the transmission prop- 
erties of a one-dimensional photonic crystal composed of 
two materials with positive and negative refractive in- 
dices. We have demonstrated unusual angular depen- 
dencies of the beam transmission through such composite 
structures which can be used for efficient beam reshap- 
ing. We have calculated the Goos-Hanchen shift of the 
refiected beam which is shown to increase dramatically 
for the angles of incidence close to the transmission band. 
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